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SUMMARY 

This report describes the development of an analytical model of a small 
turboshaft engine designed for helicopter propulsion systems. The model equa- 
tions were implemented on a hybrid computer system to provide a real-time non- 
linear simulation of the engine performance over a wide operating range. The 
real-time hybrid simulation of the engine was used as a test bed for evaluat- 
ing a microprocessor-based digital control module. This digital control module 
was developed as part of an advanced rotorcraft control program that has the 
long-term goal of rotorcraft flight and propulsion control integration. After 
tests with the hybrid engine simulation the digital control module was used to 
control a real engine in an experimental program. 

A hybrid simulation of the engine's electrical-hydromechanical control 
system was also developed. This provided a convenient way of varying the fuel 
flow and torque load inputs to the hybrid engine simulation for simulating 
transient operation. A comparison of steady-state data obtained from the 
simulation and the experimental tests is presented that validates the simula- 
tion results. Analytical model equations, analog computer diagrams, and a 
digital computer flow chart are given in the appendixes. 


INTRODUCTION 

New mission requirements for helicopters are leading to demands for faster 
response from the helicopter propulsion system and its associated control. The 
most common helicopter propulsion system employs one or two small turboshaft 
engines. Current turboshaft engine controls use pneumatic, hydraulic, and 
mechanical components to perform the control sensing, computation, and actua- 
tion. There are limitations on the dynamic response of these types of com- 
ponent and the complexity of the control functions that they can handle. One 
means of overcoming these limitations is to use digital electronic components 
to implement the control and logic functions in the control system. Digital 
control systems have proved to be feasible for other complex gas turbine con- 
trol applications (refs. 1 and 2). 

The design of helicopter propulsion controls is complicated by the very 
close dynamic coupling between the airframe and the propulsion system. Dynamic 
interaction between the fuel control system and the flexible rotor and drive- 
train has been minimized by making the fuel control very slow. Advanced, inte- 
grated controls are required to allow faster engine response while avoiding 
adverse coupling between engine and airframe. This will allow tighter rotor 
speed control and will reduce pilot workload. 



In developing helicopter engines and their control systems there are 
usually several iterations before the design is finalized. This build-and- 
test cycle can be very costly and time consuming, especially if major changes 
must be made to the control system hardware. Therefore, better analytical 
tools are needed by the control system designers. Computer simulations of gas 
turbine engines have proved to be a valuable aid in the control design process 
(ref. 3). In particular, real-time engine simulations are a cost-effective 
way of evaluating and validating digital control hardware and software prior 
to full-scale engine testing (ref. 4). 

A program is now under way at the NASA Lewis Research Center to develop 
and evaluate advanced control modes for small turboshaft engines. An initial 
task of this program was to design and build a microprocessor-based digital 
control module that could be programmed to perform the control and logic func- 
tions for the engine control system. This digital control module replaces the 
original control system hardware and performs all of the engine control func- 
tions except fuel metering and actuation of the fuel-metering valve and com- 
pressor variable geometry. For experimental engine tests the digital control 
module provides electrical signals to a modified hydromechanical unit (HMU) 
that performs the fuel metering and variable-geometry control. 

Initially the digital control module was programmed to duplicate the con- 
trol and logic functions of the original control system. To evaluate this 
digital control implementation, before actual engine testing, a real-time 
hybrid computer simulation of the small turboshaft engine was developed. This 
engine simulation provided realistic inputs to the digital control module over 
the normal operating range. After the digital control implementation was 
validated, the digital control module was used to control an actual engine. 
Later in the program, advanced control modes will be programmed on the digital 
control module and evaluated by using the real-time engine simulation prior to 
engine tests. 

To aid in validating the steady-state and transient performance of the 
engine simulation, a real-time hybrid computer simulation of the original con- 
trol system was developed. This provided a convenient way of varying the 
engine set point and facilitated running transients to assess engine simulation 
dynamics. This report documents the hybrid computer simulations of the engine 
and control system and the procedures used in developing them. 

In the following sections the engine and the control system are described. 
The mathematical modeling of engine component performance and how this was 
implemented on a hybrid computer are then discussed. The simulation of the 
original engine control system on the hybrid computer is also discussed. 
Finally, results are presented from comparisons of experimental and simulation 
steady-state data. The appendixes contain lists of equations for the engine 
and control simulation models. Detailed analog diagrams and a flow chart for 
the digital function generation program are also included in the appendixes. 


ENGINE DESCRIPTION 

The engine selected for this simulation is a small, lightweight turboshaft 
engine of the 1500-horsepower class (fig. 1). It has a five-axial-stage, 
single-centrifugal-stage compressor with variable inlet guide vanes and vari- 
able stator vanes for the first two stages. The combustor is annular with 


central fuel injectors and vortex air swirlers. The gas generator turbine 
that drives the compressor rotor is a two-stage, air-cooled, axial design. 

Air is bled from the compressor exit to cool this turbine. The two-stage power 
turbine is uncooled. It has a coaxial driveshaft that extends forward through 
the gas generator and connects through a splined joint to the engine output 
shaft. 


DEVELOPMENT OF ENGINE MODEL 

In developing the computer simulation of the small turboshaft engine, 
each of the major components (e.g., compressor, combustor, and turbine) was 
modeled by functions representing its steady-state performance data (maps) or 
a set of algebraic equations. The total volume associated with the flow pas- 
sages in the engine was distributed among a number of mixing volumes that con- 
nected the individual component models. Angular momentum equations represent 
the differential torque-speed relationship of the two rotating shafts. These 
shafts connect the compressor and gas generator turbine models and the power 
turbine and load (rotor) models. 

A block diagram of the small turboshaft engine is shown in figure 2 . The 
principal engine variables - pressures, temperatures, flows, torques, and 
rotational speeds - are shown on the diagram. The numbers are engine station 
designations. This block diagram can be viewed as a computational flow dia- 
gram. Each engine variable is represented by a line with a directional arrow 
to indicate the flow of information from one component model to another. The 
input-output relationships for a particular component are governed by thermo- 
dynamic or aerodynamic processes. These relationships are represented math- 
ematically by differential equations, algebraic equations, or single or 
multivariable functions. 

Most simulations have special requirements that depend on their applica- 
tion. For a computer simulation that is to be used for hardware-in-the-loop 
control system studies, one requirement is that it compute both the steady- 
state values and transient response of the principal engine variables over the 
operating range of interest. A second requirement is that it compute engine 
variables on a real-time basis. The latter requirement is to allow interfacing 
with the control system hardware, which operates in real time. 

After drawing the computational flow diagram the next step was to obtain 
steady-state values of the engine variables over the desired operating range. 
There were two sources of these data. The primary source was the manufac- 
turer's engine performance computer program, which was developed during the 
engine design phase. The second source was engine experimental test data. 

For the small turboshaft engine simulation development the performance program 
was used to obtain steady-state values of the engine variables for 60 operating 
points. Each operating point was defined by a value of power turbine speed 
(85 to 105 percent of design speed) and a value of power output (6.67 to 100 
percent of design power). All operating points that were selected corresponded 
to sea-level static operation. 

Mathematical equations (or tabular function data) to represent the input- 
output relationships for each engine component were determined from these 
steady-state data. The basic equations, their derivations, and the forms of 
the functions have been presented in many engine simulation reports (e.g., 
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refs. 5 to 7). For the small turboshaft engine some of the basic equations 
and functions were simplified. 

The equations and functions were first programmed in the Continuous System 
Modeling Program (CSMP) simulation language and run in nonreal time on the IBM 
370/3033 mainframe computer. The purpose was to determine if a simplified set 
of equations and functions could adequately represent the steady-state perform- 
ance of the engine. Several simplifications were considered, check cases were 
run, and steady-state results were compared with results from the more detailed 
performance program. A final model was selected that computed key engine vari- 
ables within 1 percent of the performance program values. 

The next step was to compare the results from the CSMP program with avail- 
able experimental test data. However, two problems were encountered. First, 
the engine used for the experimental tests was a nonstandard design containing 
modified compressor and turbine components. Thus the data used to develop the 
CSMP model did not represent the modified engine. Second, the experimental 
test data were insufficient to completely define the engine model. Only 12 
engine variables were measured. These included a total engine airflow 
"measurement" that was actually calculated from pressure and temperature meas- 
urements upstream of the engine. No measurements of pressures or temperatures 
were made in the hot section of the engine except for the interturbine temper- 
ature (station 45) and the exit pressure and temperature (station 49). 

To develop a CSMP model that could duplicate the experimental engine test 
data, a more complete set of engine data had to be constructed. A procedure 
for synthesizing missing engine data had been developed as an application of 
data analysis software written for a small digital computer (ref. 8). The 
synthesis procedure was based on energy and flow balance relationships for 
this type of engine. By using the 12 measured variables, 40 additional engine 
variables (and combinations thereof) were calculated at each of six test con- 
ditions. From this set of calculated and measured engine variables a modified 
CSMP model of the experimental engine was derived. The equations for the 
modified CSMP model are given in appendix B. Appendix B also briefly describes 
the mathematical representation of each component's performance. Symbols are 
defined in appendix A. 

The modified CSMP model, like most other digital simulation models, did 
not run in real time. However, it did provide certain advantages. For 
example, changes in the model structure and component performance functions 
were easily made with CSMP. Also, floating-point data obtained from the CSMP 
model, over a wide range of operating conditions, provided scaling information 
needed for the fixed-point hybrid computer simulation. 


HYBRID COMPUTER SIMULATION OF ENGINE 

The equations and functions for the experimental engine model, given in 
appendix B, were implemented on one of the Lewis Research Center's two hybrid 
computing systems. Each system consists of an EAI PACER 100 digital computer 
and two EAI 681 analog computers plus peripherals. For the engine simulation 
the digital computer was used for initial setup (i.e., setting coefficient 
potentiometers and function generators), static operation check (i.e., patch- 
ing and component checking), and digital function generation for the dual- 
variable compressor map and one of the single-variable turbine functions. One 


of the analog computers was used for all of the other calculations (i.e., sum- 
mation, integration, multiplication, division, and single-variable function 
generation). The second analog computer was used to implement the simulation 
of the original engine control system. The control system simulation is dis- 
cussed later. Analog-to-digital (A/D) and digital-to-analog (D/A) converters 
were used to transfer signals between the analog and digital computers. The 
update time of the digital computer (i.e., A/D conversion, function table 
lookup and interpolation, and D/A conversion) was less than 4 milliseconds. 
This allowed the simulation to run in real time. Diagrams for the analog com- 
puter part of the hybrid simulation and a flow chart for the digital program 
are contained in appendix C. 

To validate the hybrid computer simulation of the engine, the simulation 
was run to obtain both steady-state and transient data. Steady-state data 
were obtained by setting fixed values for inlet pressure P2, inlet temperature 
12 , outlet pressure P49, fuel flow WF, and load torque TORQLD and letting the 
computer calculate all other engine variables. However, to obtain transient 
data it was necessary to provide WF and TORQLD inputs as functions of time. 
Instead of scheduling these inputs as functions of time for a number of dif- 
ferent "real world" transients, a more direct approach was selected. A com- 
puter model of the actual engine control hardware was developed to provide the 
necessary time-varying inputs to the engine simulation in response to changes 
in the control system inputs. This model also provided information about 
internal variables in the control system to aid in developing and debugging 
software for the microprocessor-based digital control module. The engine con- 
trol system is described in the following section. 


CONTROL SYSTEM DESCRIPTION 

The engine control system was designed to provide complete and safe con- 
trol of the engine and to reduce the pilot's workload. Some of the functions 
it provides are automatic startup, compressor bleed and variable-geometry con- 
trol, rotor speed control for varying load conditions, and torque matching for 
multiple-engine installations. It also protects the engine from overspeed, 
overtemperature, and other unsafe operating conditions. 

The engine control system consists of an electrical control unit (ECU) 
and a hydromechanical control unit (HMU) . Basically, the HMU meters fuel flow 
to the engine as a function of gas generator speed NG, inlet temperature T2, 
compressor outlet pressure P3, an electrical trim signal from the ECU, plus 
pilot inputs of power-available spindle angle PAS, and load-demand spindle 
angle LDS. Rapid engine transient response to load-demand changes is provided 
by a feed-forward signal from the LDS input. Acceleration limiting, stall and 
flameout protection, gas generator speed limiting, and variable-geometry actu- 
ation are also provided by the HMU. 

The ECU provides an electrical trim signal to the HMU to satisfy the load 
requirements so as to maintain constant rotor speed and load sharing for 
multiple-engine operation. It also limits power turbine inlet temperature and 
total delivered torque and provides a separate signal to a fuel-sequencing 
valve for power turbine overspeed protection. 

Detailed information was obtained from the engine manufacturer to develop 
a model of the engine control system. A block diagram that represents the 
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principal components of the engine control system is shown in figure 3. The 
equations, functions, and logic relationships that mathematically represent 
the performance of these components are given in appendix B. Certain functions 
of the actual control system were not included in the simulation. These were 
automatic startup, variable-geometry control, and torque matching for multiple- 
engine installations. 


HYBRID COMPUTER SIMULATION OF CONTROL SYSTEM 

The equations, functions, and logic representing the control system model 
were implemented on the second analog computer of the hybrid computer system. 
The bivariable functions representing the acceleration schedule and the thermo- 
couple harness time constant were calculated on the digital computer. All 
other equations and logic were calculated on the analog computer. Analog 
diagrams and digital flow charts of the control system are included in 
appendix C. 


RESULTS AND DISCUSSION 

To determine how well the CSMP and hybrid simulations represent the engine 
steady-state performance, six test conditions were evaluated. At each test 
condition, values of inlet pressure P2, inlet temperature T2, outlet pressure 
P49, fuel flow WF, and load torque TORQLD were determined from experimental 
data. Table I lists the values for the six selected test conditions. These 
values were used as inputs to the simulations, which were run until steady- 
state conditions were reached. For the CSMP simulations all calculated vari- 
ables were printed out. For the hybrid simulations only the principal engine 
variables were recorded. 

The experimental test values and simulation results for key engine vari- 
ables are shown in table II. Agreement between simulation and experimental 
test values is quite good, with a maximum difference of about 3 percent. The 
largest differences occurred in the hot-section temperatures T45 and T49. 

Since the digital and hybrid engine simulations were programmed from the same 
model equations and function data, the observed differences between the digital 
and hybrid simulation results were attributed to slight inaccuracies in the 
analog computer components (i.e., function generators, multipliers, and 
summers) . 

To check the transient operation, the combined hybrid simulation of the 
engine and control system was exercised by manipulating the LDS and PAS inputs 
as would the pilot of a helicopter. Typically the LDS input varies the col- 
lective pitch of the helicopter rotor system and thus the torque load on the 
power turbine shaft. Changing the LDS input also sends to the fuel flow con- 
trol loop an anticipatory signal that indicates the load is changing. 

Rapid linear changes were made in the LDS input between idle and high 
load demand. Internal variables in the control system as well as principal 
engine variables were monitored. This provided a check on operation of the 
selection logic, which modifies the engine control according to various sched- 
ules and engine variable limits. The PAS input was also varied to check its 
effect on engine operation. 
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Although transient performance data were obtained from the hybrid simula- 
tion, no satisfactory transient data were available from experimental engine 
tests for comparison. However, the transient data from the simulation could 
be judged qualitatively from the standpoint of reasonableness and did provide 
helpful scaling information for developing the software for the microprocessor- 
based digital control module. 

Typical transient response of certain engine variables to changes in the 
IDS input are shown in the strip-chart recording in figure 4. The IDS input 
was ramped from 60° to 0° in 2 seconds. After a short settling period LDS was 
ramped back to 60°. The control maintained the power turbine speed NP to 
within 2 percent of the desired value. The torque reference, which represents 
the load on the power turbine, is a second-order function of LDS and NP . The 
scales on the speed traces are percent of design. The scales on the torque, 

WF, and PS3 traces are percent of full scale. The scale on the T45 trace is 
percent of the overtemperature limit value. 

After the digital control module was programmed, it was connected to the 
hybrid engine simulation. The LDS input was varied as it was for the tran- 
sients with the hybrid control simulation. The initial transient data indi- 
cated software errors. These errors were found and corrected. Subsequent 
transient response data, shown in figure 5, closely match those shown in 
figure 4. 

The digital control module was also programmed to provide logic for engine 
startup, engine overspeed protection, and sensor failure detection. All of 
these modes of control operation were successfully validated prior to actual 
engine testing with the aid of the hybrid engine simulation. 


CONCLUDING REMARKS 

A real-time hybrid computer simulation of a small turboshaft engine and 
its control system has been developed. For steady-state operation, engine 
performance data from the simulation compared very well with similar data from 
actual engine tests. Because of lack of experimental transient data no com- 
parison with the simulation data was possible. However, the simulation is 
believed to represent real transient performance well enough to be used as a 
design aid and test bed for developing a digital control module for a small 
turboshaft engine. The simulation has been used successfully to test new 
digital control hardware and software. Software errors have been found and 
corrected before actual engine testing. The simulation has provided valuable 
insights into the steady-state and transient behavior of the engine and its 
control and has therefore served as a training tool for the controls engineer. 
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APPENDIX A 


Bl 

B2 

DEL2 

DH 

DHQTH4 

DHQTH5 

dt 

EFFB 

FAR 

H 

HVF 

I 

ICX 

KDPB 

KINR 

KTRQ 

KV 

KW 

LDS 

NG 

NP 

P 

PCNG 

PR 

PS 

PS3Q2 

TH2 

T 

TAU 

THTA 

TORQ 


SYMBOLS 

Variables and Coefficients Used in Engine Simulation 

bleed fraction (function of PCNGC) 

bleed fraction (function of WA2C) 

ratio of total inlet pressure to sea-level pressure 

enthalpy drop, J/kg (Btu/lbm) 

enthalpy drop parameter for gas generator turbine, J/kg (Btu/lbm) 

enthalpy drop parameter for power turbine, J/kg (Btu/lbm) 

differential time, sec 

combustor efficiency, fractional percent 

fuel-air ratio 

enthalpy, J/kg (Btu/lbm) 

heating value of fuel, J/kg (Btu/lbm) 

p p 

polar moment of inertia, N-m-sec' 1 (ft-lbf-sec ) 
initial condition of state variable X 

2 2 2 4 

combustor pressure drop coefficient, N -sec /kg -cm -K 
( 1 bf ^-sec^l 1 bm^-i n 4 -°R ) 

angular velocity conversion coefficient, rpm-sec/radi an 
torque conversion coefficient, N-m-rpm-sec/J-radian 
( Ibf-ft-rpm-sec/Btu-radian) 

9 p © 

volume coefficient, N/cm-kg-K ( 1 bf / in -Ibm- R) 

p p 

flow coefficient, kg-cm /N-sec (Ibm-in /lbf-sec) 
load-demand spindle angle, deg 
rotational speed of compressor, rpm 
rotational speed of power turbine, rpm 

p 

total pressure, N/cm (psia) 

rotational speed of compressor, percent of design speed 
pressure ratio 

9 

static pressure, N/cm (psia) 

pressure ratio (compressor exit static to inlet total) 

ratio of inlet total temperature to standard-day temperature 

total temperature, K (°R) 

time constant, sec 

squared critical velocity ratio 

torque, N-m (ft-lbf) 
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T3Q2 temperature ratio (compressor exit total to inlet total) 
W mass flow rate, kg/sec (lbm/sec) 

WS stored mass, kg (lbm) 

WXQ2 bleed coefficient (function of WA2C) 

Qualifiers and Station Designations 


A air 

B combustor 

BL bleed air 

C corrected value 

CM compressor 

DSN design value 

F fuel 

GT gas turbine 

L lagged 

LD load 

PH per hour 

PT power turbine 

2 compressor inlet 

25 intercompressor bleed port 

3 compressor outlet 

31 combustor inlet 

41 gas generator turbine inlet 

44 gas generator turbine outlet 

45 power turbine inlet 

49 power turbine outlet 

Variables and Coefficients Used in Control System Simulation 

DWFPL incremental fuel flow demand signal due to LDS 

ICX initial condition of state variable X 

KGVNR NP governor rate compensation gain 

KINTG integral gain (ECU compensation) 

KNDRP NG droop line slope 

KNGVN additional NP governor gain 

KNTRM speed error trim gain (normal), V/rpm 
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KPROP 

KSPDS 

KTMC 

KTMFB 

KTMLO 

KTMRU 

KTMSR 

KTSG1 

KT45LM 

KWFD 

KWLDSR 

PAS 

PCNGD 

PCNGDH 

PCNGLDL 

PCNGI 

PCNGS 

PCNPRF 

PCNPS 

PS3L 

RULDS 

SPDB 

SPDER 

SPDG 

SPDL 

SPDR 

SPRDC 

SPRD2 

SPDS 

SPDSF 

SPDSHL 

SPDSI 

SPDS IL 

SPDSL 

SPDSLL 


proportional gain (ECU compensation) 

NP governor gain 

torque motor coil conductance, mA/V 

torque motor LVDT feedback gain, V/cm (V/in.) 

torque motor sensitivity, cm/sec-mA (in./sec-mA) 

torque motor linkage gain, cm/cm (in. /in.) 

torque motor compensation gain, V/V 

T45 limiter compensation gain 

T45 limiter gain 

midrange slope (deceleration schedule) 
reference level (LDS WF/PS3 demand signal) 
power-available spindle angle, deg 
NG demand based on LDS, percent of design speed 
upper NG limit (deceleration schedule), percent of design speed 
lower NG limit (deceleration schedule), percent of design speed 
reference for idle NG droop line (function of T2), percent of design 
speed 

sensed gas generator turbine speed, percent of design speed 
power turbine speed reference, percent of design speed 
sensed power turbine speed, percent of design speed 
sensed PS3, N/cm (psia) 

fuel flow demand trim signal from LDS and ECU 

speed error threshold for additional NP governor gain, percent of 
design speed 

power turbine speed error, percent of design speed 
ECU trim demand signal, V 

speed error with governor rate lag compensation 

speed error above threshold value, percent of design speed 

speed error with rate compensation 

speed error with governor rate compensation 

speed error with governor lag 

speed error with governor dynamics 

high limit of integral of speed error 

speed error with integral compensation 

integral of speed error (after limiter) 

speed error with governor lead-lag 

low limit of integral of speed error 
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SPDSP speed error with proportional compensation 
SPDSS selected maximum error (NP or T45 overlimit) 

SPDS1 speed error with trim gain and governor rate compensation 

SPDS2 speed error with additional gain and governor rate compensation 

TCGRLD time constant (NP governor rate compensation, lead), sec 

TCGRLG time constant (NP governor rate compensation, lag), sec 

TCGVLD time constant (governor loop lead), sec 

TCGVLG time constant (governor loop lag), sec 

TCNGS time constant (NG sensor lag), sec 

TCPS3 time constant (PS3 sensor), sec 

TCSPDG time constant (ECU trim signal lag), sec 

TCTMSR time constant (torque motor lag compensation), sec 

TC45 time constant (T45 harness lag (fixed)), sec 

TC45H time constant (T45 harness lag (variable)), sec 

TC45L time constant (T45 limiter lag), sec 

TC45LD time constant (T45 limiter lead compensation), sec 

TC45LG time constant (T45 limiter lag compensation), sec 

TERR T45 overtemperature error, K (°R) 

TMBS1 torque motor error bias, V 

TMBS2 torque motor current bias, mA 

TMDB torque motor current deadband, mA 

TMI torque motor current, mA 

TMILIM torque motor current limit, mA 

TMLO torque motor flapper position, cm (in.) 

TMLDL torque motor flapper position limited, cm (in.) 

TMRU WF/PS3 trim signal 

TMSE torque motor error signal, V 

TMSIG torque motor current (with deadband), mA 

TMSL intermediate torque motor error compensated, V 

TMSLL torque motor error compensated, V 

TMSR torque motor error compensated, V 

TSGE T45 error with lag, K (°R) 

TSGL intermediate T45 error compensated, K (°R) 

TSIG T45 error compensated, K (°R) 

TSIGS T45 positive error signal to selector, K (°R) 

T45C0R T45 harness correlation, K (°R) 

T45E T45 plus harness correlation, K (°R) 
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T45EL T45 sensed, K (°R) 

T45L T45 sensed intermediate, K (°R) 

T45LIM T45 maximum limit, K (°R) 

WF fuel flow (after limiting), kg/hr (lbm/hr) 

WFIDM WF/PS3 demand limit signal (idle schedule) 

WFIRF WF/PS3 demand value (idle schedule limiting value; function of T2) 

WFMAX maximum fuel flow limit, kg/hr (lbm/hr) 

WFMIN minimum fuel flow limit, kg/hr (lbm/hr) 

WFMV fuel flow (after limiting), kg/hr (lbm/hr) 

WFPAC WF/PS3 demand limit signal (acceleration schedule) 

WFPDC WF/PS3 demand limit signal (deceleration schedule) 

WFPDCH upper limit WF/PS3 demand (deceleration schedule) 

WFPDCL lower limit WF/PS3 demand (deceleration schedule) 

WFPDM combined WF/PS3 demand signal before limits 
WFPDRP AWF/PS3 demand signal (droop line) 

WFPDO WF/PS3 demand signal (after idle limit selector) 

WFPD1 WF/PS3 demand signal (after acceleration limit selector) 

WFPD2 WF/PS3 demand signal (after acceleration limit selector) 

WFPRF AWF/PS3 demand signal (function of PAS) 

WFPTP WF/PS3 demand value (topping level) 

WFQD AWF/PS3 demand signal (function of LDS) 

W45R power turbine flow parameter, kg-cm^-K^ ^/N-sec (lbm-in^-°R^/lbf-sec) 
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APPENDIX B 


MATHEMATICAL REPRESENTATION OF ENGINE AND CONTROL COMPONENTS 

Engine Components 

Compressor. - The corrected compressor airflow is calculated from an 
overall performance map that is a bivariable function of pressure ratio and 
corrected rotor speed. The temperature ratio across the compressor is a non- 
linear function of pressure ratio. The bleed flow fractions, for compressor 
interstage bleed, are functions of corrected rotor speed and airflow. 


PCNGC = NGC/NGDSN * 100 

(Bl) 

PS3Q2 = PS3/P2 

(B2 ) 

WA2C = f 2 (PS3Q2, PCNGC) 

(B3) 

WA2 = WA2C * DEL2//TH2 

(B4) 

H2 = 889.5 * T2 

(B5) 

T3Q2 = f 2 (PS3Q2) 

(B6) 

T3 = T3Q2 * T2 

( B7 ) 

B1 = f 3 (PCNGC) 

(B8) 

B2 = f (WA2C) 

(B9) 

WBL25 = ( Bl + B2 ) * WA2 

(BIO) 

WA3 = WA2 - WBL25 

(Bll) 


Volume 3 . - Two dynamic equations and a state equation are used to cal- 
culate the total pressure, temperature, and stored mass (air) in volume 3. 

The bleed fraction is a function of corrected compressor airflow. A first- 
order lag equation is used to calculate temperature in this volume. Enthalpy 


is calculated as a linear function of temperature. 

WXQ2 = f 5 (WA2C) (B12 ) 

WBL3 = (WXQ2 + 0.0025) * WA2 (B13) 

WS3 = J (WA3 - WBL3 - WA31)dt + ICWS3 (B14) 

o 

ft 

T3L = J 1/TAU3 * ( T3 - T3L)dt + ICT3 (B15) 

o 
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H3 = 928.9 * T3L - 16538.5 
P3 = KV3 * WS3 * T3L 


(B16) 
( B 1 7 ) 


Combustor. - The combustor mass flow is a function of inlet pressure, 
inlet temperature, outlet pressure, and a pressure loss coefficient. The com- 
bustor efficiency is a function of fuel-air ratio. The enthalpy increase 
through the combustor is a function of fuel-air ratio, combustor efficiency, 
and fuel properties. 


WA31 = [P3 *(P3 - P41)/(KDPB * T3)] 1/2 (B18) 

WF = WFPH/3600 (B19) 

FAR = WF/WA31 (B20) 

EFFB = f g (FAR) (B21) 

DHB « FAR * EFFB * HVF (B22) 

H41 = ( H3 + DHB) / ( 1 . + FAR) (B23) 


Volume 41 (combustor outlet). - The temperature in volume 41 is calcu- 
1 ated as a linear function of enthalpy. The state equation is combined with a 
dynamic continuity equation to calculate pressure. 

T 41 = 8.926E-4 * H41 + 180.4 (B24) 


p 41 = J K V 41 * T 41 * (WA31 + WF - W41)dt + ICP41 (B25) 

o 

Gas generator turbine . - In prior simulations (refs. 5 to 7) turbine per- 
formance was represented by an enthalpy drop parameter and a flow parameter 
that were both bivariable functions of pressure ratio and corrected speed. 
However, for this engine the turbine data obtained from the manufacturer s 
engine performance program were worked into a simpler functional form. By 
using a critical velocity parameter (ref. 9), an enthalpy drop parameter was 
obtained that was a function of only pressure ratio. From these parameters 
the actual enthalpy drop, exit enthalpy, and temperature can be calculated. 
Because over the operating range of interest this turbine was found to be 
choked, a simple choked nozzle equation was used to calculate flow. 


PRGT = P45/P41 (B26) 

DHQTH4 = f 7 (PRGT) (B27) 

THTA41 = 2.9322E-3 * T41 + 0.0856 (B28) 

DHGT = DHQTH4 * THTA41 (B29) 

H44 = H41 - DHGT (B30) 


!ii i r 


14 



W 41 = KWGT * P41/THTA41 


(B31 ) 


Volume 45 (interturbine) . - In volume 45 bleed cooling flow is mixed with 
gas generator turbine exit flow. The enthalpy of the mixed gases is assumed 
to be proportional to the enthalpy of the turbine exit gas. The temperature 
is calculated as a linear function of enthalpy. The state equation is combined 
with a dynamic continuity equation to calculate pressure. 

H45 = 0.96 * H44 (B32) 

T45 = 9.456E-4 * H45 + 111.9 (B33) 

P 45 = J KV45 * T 45 * (W41 - W45 + 0.7826 * WXQ2 * WA2) (B34) 

o 

Power turbine . - For this turbine an enthalpy drop parameter and a flow 
parameter are obtained that are functions of only the pressure ratio across 
the turbine. By using these parameters and the critical velocity parameter 
for this turbine, the enthalpy drop, exit enthalpy, exit temperature, and 
weight flow can be calculated 


PRPT = P49/P45 

( B35 ) 

DHQTH5 = fg (PRPT) 

(B36 ) 

THTA45 = 2.9322E-3 * T45 + 0.0856 

(B37 ) 

DHPT = DHQTH5 * THTA45 

( B38 ) 

H49 = H45 - DHPT 

(B39 ) 

T 49 = 9.448E-4 * H49 + 107.7 

(B40) 

W45C = f 9 (PRPT) 

(B41 ) 

W45 = W45C * P45/THTA45 

(B42 ) 


Compressor and gas generator turbine rotor . - The torques of these com- 
ponents are calculated from the changes in the fluid energy from input to out- 
put. For the compressor an empirical relationship is used to account for 
interstage compressor bleed. The rotor speed is determined from a dynamic 


angular momentum equation. 

TORQCM = KTRQ/NG * [0.29 * H3 * WA3 + (0.71 * H3 - H2) * WA2] (B43) 
TORQGT = KTRQ/NG * DHGT * W41 (B44) 

t 

NG = f KINR/IGT * (TORQGT - T0RQCM)dt + ICNG (B45) 

o 

NGC = NG/ /rH2 (B46) 
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Power turbine and load rotor . - The power turbine torque and rotor speed 
are calculated by using the same types of equations as above. The load torque 
is determined as a function of load-demand spindle angle LDS and is speed 
dependent. 


TORQPT = KTRQ/NP * DHPT * W45 (B47) 

TORQLD = [50.843 - LDS * (0.0835 - 0.1018 * LDS)] * (NP/20000.) ** 2 (B48) 

_t 

NP = f KINR/IPT * (TORQPT - T0RQLD)dt + ICNP (B49) 

o 

Control Components 
Electrical Control Unit 

T45 thermocouple harness . - The dynamics of the T45 thermocouple harness 
are represented by two first-order lags. One of the time constants is fixed 
and the second is a bivariable function of the temperature and flow parameter 
at station 45. 


T45E = T 45 + T45C0R (B50) 

W45R = W45 * T45/P45 (B51) 

TC45H = f 1Q (W45R, T45) (B52) 

t 

T45L = / 1/TC45 * (T45E - T45L)dt + ICT45E (B53) 

o 

t 

T45EL = / 1/TC45H * (T45L - T45EL)dt + ICT45E (B54) 

o 


T 45 limit control. - In this control the sensed temperature T45 is com- 
pared with a limit value to give an error signal. This signal then goes 
through a first-order lag and a lead-lag circuit to give a compensated T45 
error signal . 


TERR = T45EL - T45LIM 


(B55 ) 

t 

TSGE = f KT45LM/TC45L * (TERR - TSGE)dt + 
0 

ICTSGE 

( B56) 

t 

TSGL = / KTSG1/TC45LG * (TSGE - TSGL)dt + 
0 

ICTSG1 

(B57) 

TSIG = TC45LD * (TSGE - TSGL)/TC45LG 


(B58) 


NP speed control. -In this control the sensed speed NP is compared with a 
reference value to given an error signal. This signal then goes through a 


■ in r 
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lead-lag and first-order lag to give a compensated NP error signal. Additional 
gain and error rate compensation is added if the error signal is larger than a 
small threshold value. 


SPDER = PCNPS - PCNPRF (B59) 

(SPDER - SPDB for SPDER > SPDB 

SPDR = <0.0 for -SPDB < SPDER < SPDB (B60) 

(SPDER + SPDB for SPDER < -SPDB 

t 

SPDL = f 1/TC6RLG * (SPDR - SPDL)dt + 0.0 (B61) 

o 


t 

SPDRC = f 1/TCGRLD * (SPDL - SPDRC) dt + 0.0 (B62) 

o 

SPDR2 = KGVNR * (SPDL - SPDRC) /TCGRLD (B63) 

SPDS1 = KNTRM * SPDER + KNGVN * SPDR + SPDR 2 (B64) 

SPDS2 = 2.0 * SPDS1 (B65) 


t 

SPDS = f KSPDS/TCGVLG * (SPDS2 - SPDS)dt + 0.0 (B66) 

o 


SPDSL = TCGVLD * ( SPDS2 - SPDS)/TCGVLG 


( B67 ) 


t 

SPDSF = / 1/TCGVLG * (SPDSL - SPDSF)dt + 0.0 (B68) 

o 


Speed-temperature limit selector . - This selector logic normally selects 
the speed error signal. If the temperature error indicates T45 over- 
temperature, the logic selects the error signal that will give lowest fuel 
flow. 


TSIGS = 


(tsig 

JO .0 


for TSIG > 0.0 
for TSIG < 0.0 


(B69 ) 


SPDSS = 


(tsigs 

1SPDSF 


for TSIGS > SPDSF 
for TSIGS < SPDSF 


( B70) 


Principal ECU compensation . - The ECU compensation circuit contains 
parallel proportional and integral paths followed by a first-order lag. The 
integrated error signal is limited if it exceeds maximum or minimum values. 


SPDSP = KPROP * SPDSS 


(B71 ) 


SPDSI = / (KINTG * SPDSS)dt + ICSPDI (B72) 

o 
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ISPDSHL 
SPDSIL = <SPDSI 
(SPDSLL 


for SPDSI > SPDSHL 

for SPDSLL < SPDSI < SPDSHL 

for SPDSI < SPDSLL 


(B73) 


SPDG = / 1/TCSPDG * (SPDSIL + SPDSP - SPDG)dt + 1CSPDG 
o 


(B74) 


Hydromechanical Control Unit 

The primary function of the hydromechanical control unit is to control 
the fuel flow to the engine by positioning the fuel-metering valve. The valve 
position is affected by a combination of pilot-initiated inputs, sensed engine 
variables, and an electrical trim signal from the ECU. The valve position is 
also constrained by an idle schedule, an acceleration schedule, a deceleration 
schedule, and maximum and minimum absolute fuel flow limits. In the following 
description the fuel flow demand parameter WFPDM is the primary control vari- 
able. Fuel valve position and fuel flow are proportional to this parameter 
multiplied by the sensed compressor discharge pressure PS3L. 

Topping line. - A maximum value of the fuel flow demand parameter is cal- 
culated as a function of compressor inlet temperature T2. 

WFPTP = f n (T2 ) (B75) 

Gas generator turbine speed and droop line . - An incremental fuel flow 
demand parameter is calculated from the difference between sensed speed NG and 
101 percent of design speed. This provides NG droop compensation through pro- 
portional feedback control. 


PCNG = NG/NGDSN * 100 (B76) 

t 

PCNGS = f 1/TCNGS * (PCNG - PCNGS)dt + ICPCNG (B77) 

o 

WFPDRP = KNDRP * (101. - PCNGS) (B78) 

Pow er-available spindle input . - This pilot-initiated input determines 
another incremental fuel flow demand parameter. This in effect sets the maxi- 
mum available gas generator turbine speed and the maximum available power. 

WFPRF = f 12 (PAS) ( B79 ) 

Load-demand spindle input. - This pilot-initiated input changes the col- 
lective pitch of the aircraft rotor system and thus changes the load on the 
power turbine. It also provides a feedforward signal (effectively another 
incremental fuel flow demand parameter) to rapidly change NG in anticipation 
of the load change. 


PC.NGD = f 
WFQD = f 


13 (LDS) 

( B80) 

1 4 (LDS) 

(B81 ) 
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DWFPL = KWLDSR - WFQI + KNDRP * (101. - PCNGD) 


( B82 ) 


Torque motor. - The torque motor receives a trim demand signal from the 
ECU. This signal is compared with a position feedback signal to produce an 
error signal. The error signal goes through a first-order lead-lag circuit 
and is amplified to provide a voltage that is applied to the torque motor coil. 
The resulting torque motor current moves the flapper valve to a position that 
will provide a null current (within a deadband) through the position feedback 
signal. The flapper valve movement is constrained by position limits and is 
amplified mechanically to provide a fuel flow demand trim signal. 


TMSE = SPDG - KTMFB * TMLO - TMBS1 (B83) 

J. 

TMSR = f KTMSR/TCTMSR * (TMSE - TMSR)dt + ICTMSR (B84) 

o 

TMSL = TMSE + TMSR (B85) 

TMSLL + KTMC * TMSL (B86) 

(TMILIM for TMSLL > TMILIM 

TMI = < TMSLL - TMBS2 for -TMILIM < TMSLL < TMILIM (B87) 

(-TMILIM for TMSLL < -TMILIM 

(TMI - TMDB for TMI > TMDB 

TMSIG = (0.0 for -TMDB < TMI < TMDB (B88) 

(TMI + TMDB for TMI < - TMDB 

TMLO = f (KTMLO * TMSIG)dt + ICTMLO (B89) 

o 

(TMLMAX for TMLO > TMLMAX 

TMLOL = (TMLO for TMLMIN < TMLO < TMLMAX (B90) 

(TMLMIN for TMLO < TMLMIN 

TMRU = KTMRU * TMLOL (B91) 


Fuel flow demand parameter. - A combined fuel flow demand parameter is 
obtained by summing the topping line value and the incremental values from the 
other inputs. The feedforward signal from the LDS input and the trim signal 
from the torque motor are summed, and the sum is limited to positive values. 
Because this signal is subtracted from the other inputs, it can only trim down. 


RULDS = 


(dWFPL + TMRU 

( 0.0 


for (DWFPL + TMRU) >0.0 
for (DWFPL + TMRU) <0.0 


(B92) 


WFP DM = WFPTP + WFPDRP - WFPRF - RULDS (B93) 


Idle schedule . - This schedule determines a limiting value of fuel flow 
demand parameter for PAS in the ground-idle region (23.5 < PAS < 28.5). This 
parameter is a function of NG and T2. 
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(B94) 


WFIRF = f 15 (T2 ) 

PCNGI = f lg (T2) ' (895 ) 

WFIDM = KNDRP * (PCNGI - PCNGS) + WFIRF (B96) 

Acceleration schedule. - This schedule determines a limiting value of 
fuel flow demand parameter during engine acceleration. This parameter is a 
bivariable function of NG and T2. 

WFPAC = f 17 (PCNGS, T2) (B97) 

Deceleration schedule. - This schedule determines a limiting value of 
fuel flow demand parameter during engine deceleration. This parameter is a 
function of NG. 


( WFPDCH for PCNGS > PCNGDH 

WFP DC = <KWFD * (PCNGS - PCNGDL ) for PCNGDL < PCNGS < PCNGDH (B98) 

(WFPDCL for PCNGS < PCNGDL 

Select logic - fuel flow demand parameter. - This logic selects the final 
fuel flow demand parameter by considering the limiting values determined by 
the schedules. 


WFPDO = 


WFP DM 
WFIDM 


for WFPDM > WFIDM 
for WFPDM < WFIDM 


(B99 ) 


WFPD1 - 


(wFPDO 

WFPAC 


for WFPDO < WFPAC 
for WFPDO > WFPAC 


(B100) 


WFPD2 = 


(wFPDl 

(WFPDC 


for WFPD1 > WFPDC 
for WFPD1 < WFPDC 


(B101 ) 


Compressor static discharge pressure sensor , - The response of this sensor 
is considered to be a first-order lag. 


t 

PS3L = f 1/TCPS3 * (PS3 - PS3L)dt + ICPS3 (B102) 

o 

Fuel-flow-metering valve. - Mechanical stops on this valve limit the fuel 
flow to maximum and minimum values. 

WFMV = WFPD2 * PS3L (B103) 


l WFMAX 
WFPH = < WFMV 
(WFMIN 


for WFMV > WFMAX 

for WFMIN < WFMV < WFMAX 

for WFMV < WFMIN 


( B104) 


APPENDIX C 


ANALOG DIAGRAMS AND DIGITAL FLOW CHART 


Flow Chart of Function Generation Program 
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TABLE I . - SFLECTEO TEST CONDITIONS 


Test 

condition 

Compressor 
inlet total 
pressure, 
P2 

Compressor 
inlet total 
temperature, 
T? 

Power tur- 
bine inlet 
pressure, 
P49 

Load torque, 
TORQLD 

Fuel flow, 
WF 

N-m 

ft-lbf 

kg/hr 

Ib/hr 

N /C(T? 

psia 

K 

°R 

N/cm ? 

psia 

1 

9.91 

14.37 

287.1 

516.7 

9.91 

14.37 

40.8 

30.1 

63.5 

140.1 

2 

9.77 

14.17 

286.4 

515.6 

9.95 

14.43 

122.2 

90.1 

134.8 

297.2 

3 

9.76 

14.16 

282.4 

508.3 

9.97 

14.46 

201.1 

148.3 

168.7 

372.0 

4 

9.71 

14.09 

282.2 

508.0 

10.07 

14.60 

280.0 

206.5 

207.9 

458.4 

5 

9.67 

14.02 

281.8 

507.2 

10.09 

14.63 

371.9 

274.3 

254.3 

560.6 

6 

9.60 

13.92 

281.8 

507.2 

10.15 

14.72 j 

489.2 

360.8 

315.0 

694.4 


TABLE II. - COMPARISON OF STEADY-STATE RESULTS 


Test 

condition 

Type of 
data 

Compressor 
outlet total 
pressure, 

P3 

Compressor 
outlet total 
temperature, 
T3 

Compressor 
inlet air 
mass flow, 
WA2 

Rotational 
speed of 
compressor , 
NG, 
rpm 

Rotational 
speed of 
power 
turbine, 
NP, 
rpm 

Power turbine 
inlet total 
temperature , 

T 45 

Power turbine 
outlet total 
temperature , 

T 49 

N/cm 2 

psia 

K 

°R 

kg/sec 

1 bm/sec 

K 

"R 

K 

°R 

1 

Experimental 

40.0 

58.0 

462 

832 

1.45 

3.20 

29 480 

10 995 

785 

1413 

747 

1345 


Digital 

39.4 

57.2 

461 

830 

1.47 

3.24 

29 481 

10 995 

773 

1392 

738 

1328 


Hybr i d 

39.6 

57.4 

461 

830 

1.47 

3.25 

29 580 

11 115 

767 

1380 

731 

1316 

2 

Experimental 

78.0 

113.1 

570 

1026 

2.34 

5.16 

37 875 

19 995 

876 

1577 

779 

1402 


Digital 

75.5 

109.5 

563 

1013 

2.36 

5.21 

37 862 

19 997 

880 

1584 

778 

1400 


Hybrid 

75.8 

110.0 

565 

1017 

2.39 

5.27 

37 955 

20 235 

868 

1563 

764 

1376 

3 

Experimental 

95.8 

139.0 

601 

1081 

2.79 

6.16 

39 236 

19 995 

903 

1626 

783 

1410 


Digital 

95.8 

139.0 

599 

1079 

2.80 

6.17 

39 239 

19 921 

900 

1620 

768 

1383 


Hybrid 

96.7 

140.3 

601 

1082 

2.83 

6.24 

39 335 

20 212 

894 

1609 

761 

1370 

4 

Experimental 

111.1 

161.1 

626 

1127 

3.14 

6.92 

40 409 

19 995 

962 

1731 

809 

1457 


Digital 

110.7 

160.6 

624 

1124 

3.13 

6.90 

40 398 

19 788 

957 

1723 

796 

1432 


Hybrid 

111.4 

161.5 

626 

1126 

3.15 

6.95 

40 465 

20 040 

956 

1721 

793 

1427 

5 

Experimental 

127.4 

184.8 

652 

1173 

3.47 

7.66 

41 393 

19 995 

1021 

1838 

849 

1529 


Digital 

126.7 

183.8 

652 

1173 

3.47 

7.65 

41 382 

19 952 

1021 

1838 

827 

1488 


Hybrid 

127.0 

184.2 

652 

1173 

3,48 

7.68 

41 405 

19 977 

1018 

1832 

823 ! 

1482 

6 

Experimental 

146.1 1 

211.9 

682 

1228 

3.86 

8.50 

42 880 

19 995 

1097 

1974 

889 

1600 


Digital 

146.1 

211.9 

683 

1229 

3,86 

8.50 

42 861 

19 998 

1099 

1979 

870 

1566 


Hybrid 

146.5 

212.5 

683 

1230 

3,87 

8.54 

42 970 

20 107 

1093 

1968 

864 : 

1556 
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Figure 1. - Small turboshaft engine. 



Figure 2. - Block diagram of small turboshaft engine. 














Figure 3, - Block diagram of engine control system. 
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